Dedicated to John T.Lewis 



Bose-Einstein Condensation in the Luttinger-Sy Model 



^ : Olivier LENOBLE^ and Valentin A. ZAGREBNOV^ 

O 

j__^ . University of California at Irvine, Department of Mathematics, 

Ph! 92697-38 Irvine, California, USA^ 

< 

Universite de la Mediterranee and Centre de Physique Theorique, 

^ ■ Luminy-Case 907, 

13288 Marseille, Cedex 09, France^ 

> 

00 

o 

O ' Abstract 
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j^ , averaged one-body reduced density matrix. We show that mathematical mechanism of 

the Bose-Einstein condensation in this random model is similar to condensation in a 
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model we prove the Kac-Luttinger conjecture, i.e., that this model manifests a type I BEC 
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1 Introduction 

1.1 In our recent paper ^2]; we presented some general mathematical results concerning the 
existence of the Bose-Einstein condensation (BEC) of the Perfect Bose-Gas (PBG) placed in a 
semi-bounded from below homogeneous ergodic random external potential (random impurities). 
There we show that for the infinite- volume one-particle Schrodinger operator, a generic Lifshitz 
tail behaviour of density of states near the lower edge of the spectrum reduces the critical 
dimensionality of the BEC for PBG from dimensionality d = 2 + etod = l. Therefore, the 
randomness enhances the BEC, and moreover, it is shown to be stable with respect to the 
mean-field particle interaction |12j. 

To tackle the corresponding Off-Diagonal Long-Range Order (ODLRO), we introduced in 
[T^ a concept of the space- averaged one-body reduced density matrix. In spite of a rather 
accurate estimate of this matrix out of the BEC domain, that shows an enhancement of the 
ODLRO exponential decay due to impurities, we did not obtain in [TH] any sound estimate for 
this order in the BEC domain. 

1.2 The aim of this paper is to present a rigorous study of a particular case of a one- dimensional 
PBG model in homogeneous ergodic non-negative random potential induced by the Poisson 
distributed singular point impurities (the Luttinger-Sy model f?! ^]). We show that this 
model allows a rigorous mathematical approach to condensation and that one can compute 
explicitly some of thermodynamical quantities even in the BEC domain. This concerns in 
particular the ODLRO behaviour of the two-point correlation function (space-averaged one- 
body reduced density matrix) in the condensation regime. 

Notice that the first study of possible modifications of d = 3 dimensional BEC in the PBG 
caused by repulsive finite-range impurities goes back to Kac and Luttinger [SI El- They predicted 
an enhancement of d = 3 BEC by indication that due to impurities there is decreasing of the 
critical density , but did not discuss a modification of the critical dimensionality. They also 
mentioned a puzzling question about the nature of the established BEC. For example, they 
conjectured that this condensate occurs as a macroscopic occupation of only the ground-state: 
type I BEC. We prove this conjecture in the case of the Luttinger-Sy model, see discussion in 
Section 6. We show that the nature of BEC in this model is close to what is known as the 
" Bose-glass" , since it may be localized by the random potential. This is of interest for example 
in experiments with liquid ^He in random environments like Aerogel and Vycor glasses, |milUj. 

On the other hand, the nature and behaviour of the lattice BEC may be quite different. 
First of all, the lattice Laplacian and the on-site Bose-Hubbard particle repulsion produces a 
coexistence of the BEC {superfluidity) and the Mott insulating phase as well as domains of 
incompressibility, see e.g. a very complete review [22] • Adding disorder makes the correspond- 
ing models much more complicated. The physical arguments show that the randomness may 
suppress the BEC (superfluidity) as well as the Mott phase in favour of the localized Bose- 
glass phase, but this is very sensitive to the choice of the random distribution, for some recent 
rigorous results see [5]. 

1.3 The paper is organized as follows. In Section 2 we recall definition of the Luttinger-Sy model 
and some of its properties. We prove the self- averaging of corresponding integrated density of 
states in Section 3 and we calculate it explicitly. In Section 4 we prove that the established 
integrated density of states implies the existence of generalized BEC in the case of PBG. 

Our main results are collected in Sections 5 and 6. There we recall the notion of the space- 
averaged one-body reduced density matrix and we prove that it has an almost sure nonrandom 
thermodynamic limit (self-averaging), which can be calculated explicitly for all values of particle 



density. We also prove that randomness enhances decay of the two-point correlation function. 
In particular we show that it keeps this decay always exponential, even in the presence of 
BEG. We found that in the latter case the ODLRO is non-zero and that it coincides with the 
condensate density. 

The properties of the BEG are discussed in concluding Section 6. First, we analyze the 
critical density dependence on the amplitude of the repulsive Poisson point impurities. Notice 
that for the Luttinger-Sy model the singular point impurities mean that this amplitude is 
infinite. Next, we study the problem of the condensate nature and its localization. To elucidate 
this point we invented a hierarchical one-dimensional nonrandom model, which mimics in a 
certain sense the (random) Luttinger-Sy model. We show that this hierarchical model can 
manifest different types (I, II and III) of generalized van den Berg-Lewis-Pule condensations 
PP localized in one, several or infinite number of {infinite) intervals of logarithmic sizes. To 
discriminate between these options, i.e. to prove or disprove the Kac-Luttinger conjecture, one 
has to have a quite detailed information about the energy level spacing in random intervals 
generated by the Poisson impurity positions. We prove this conjecture, i.e., that type I BEG in 
the Luttinger-Sy model is localized in a single "largest" (i.e. infinite) interval of the logarithmic 
size. 

2 The Luttinger-Sy Model 

In the framework of general setting this model corresponds to the following one- dimensional 
{d = 1) single-particle random Schrodinger operator in the Hilbert space 7i = L^(M): 

2.1 Gonsider a random (measurable) potential v^'\-) : f] x M — > M, {uj,x) h^ f'^(x), which is a 
random field on a probability space (f2,jF, P), with the properties: 

(a) v'^ is homogeneous and ergodic with respect to the group {r^}^^^ of probability-preserving 
translations on {Q,J-',F); 

(b) v'^ is non-negative and inf^jgigd {v'^{x)} = 0. 

By E {■} := f^ F{duj) {■} we denote the expectation with respect to the probability measure in 
{Q, T , P). Then the random Schrodinger operator corresponding to the potential v^ is a family 

of random operators {/i^^j^g^ • 

/i^:=t + t,-, (2.1) 

where t := (—A/2) is the jree one-particle Hamiltonian, i.e., a unique self-adjoint extension of 
the operator: —A/2, with domain in L^(M). 

Notice that assumptions (a) and (b) guarantee that there exists a subset f2o C JF with 
P(r2o) = 1 such that operator (J2.H) is essentially self-adjoint on domain C^(M) for every cj G fig 
(see e.g. [201 Gh.1.2). 

2.2 Let u{x) > 0, a; e M, be continuous function with a compact support. We call it a 
(repulsive) single-impurity potential. Let {^xidx)}^^^ be random Poisson measure on M with 
intensity A > : 

P ({cu e fi : z/^A) = n}) = ih}MLe-^\M ^ n G No = N U {0} , (2.2) 

n! 

for any bounded Borel set A C M. Then the non-negative random potential v'^ generated by 
the Poisson distributed local impurities has realizations 
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I vt{dy)u{x-y)= J2 ui^-x"^)- (2.3) 



Here the random set X'^ corresponds to impurity positions X^ = {x*^} . C M, which are the 
atoms of the random Poisson measure, i.e., card {X'^ \ A} = i^^(A) equals to the number of 
impurities in the set A. Since the expectation E(i^^(A)) = A |A|, the parameter A coincides 
with impurities concentration on the axe M. 

Remark 2.1 The random potential \2.'J\f is obviously homogeneous and ergodic (even strongly 
mixing), i.e. it verifies the conditions (a) and (h). Moreover, see \20^ Ch.II.5, we have that: 

- There exists a nonrandom measure dM{E) on M such that 

dM{E):=¥.{£h^{dE-Q,Q)}. (2.4) 

Here Sh'^{dE;x,y) is the kernel of the spectral decomposition measure corresponding to random 
Schrodinger operator h'^ . The spectrum o'{h'^) of h'^ is almost-surely (a.s.) nonrandom and it 
coincides with the support of M : (y{h^) = suppA/". 

- For repulsive impurities with compact support and for Poisson distribution, the a.s nonrandom 
spectrum a{h'^) = M+. Thus the lower edge of the spectrum ini {a{h'^)} = 0, i.e. it coincides 
with the lower edge of the spectrum of the nonrandom operator t, see \2.1]] . 

- In one- dimensional case the asymptotic behaviour of the integrated density of states J\f{E) : = 
A/'((— oo, -E]) as E I has the form (the Lifshitz tail): 

\nU{E)^-\{^-^y\ EiO, (2.5) 

for Q > 0. Recall that in the nonrandom case v^ = Q one obtains: N'iE) ~ E'^^'^, E I 0. 

2.3 Luttinger and Sy defined their d = 1 model ^Tj restricting the single-impurity potential 
to the point ^-potential with amplitude a > 0. In fact this choice (even for more general case 
of random {oj}) goes back to Frish and Lloyd [7j. Then the corresponding random potential 
(j2.3p takes the form: 

<{x) := I ut{dy)a5{x - y) = a J] 5{x - xj). (2.6) 

Now the self-adjoint one-particle random Schrodinger operator 

/i-:=t + <, (2.7) 

can be defined in the sense of the sum of quadratic forms. In spite of a singular nature of 
this random potential, by standard limiting arguments [20] it inherits the properties quoted in 
Remark 12.11 

2.4 Moreover, the same arguments [201 are applied to define a strong resolvent (s.r.) limit 
of Hamiltonians ()2.7|) . when a — > +cx), which is the last step in definition of the Luttinger- 
Sy model, ^7]. This limit gives the self-adjoint (Friedrichs) extension of symmetric operator 
to = -A/2 with domain dom{to) = {f e H : f e C^{R \ X'^)}. For any cu G 1^ we denote this 
extension by 

h^^ := s.r. hm h^, . (2.8) 

Since for any uj E Q the set X'^ can be ordered: X^ = {x^} . , it generates a set of intervals 
{/^ := {x'^-i, x'j)} ■ of lengths {L^ := x^ — x"^^^^ . . Then one can decompose the Hilbert space 
7i = L^(]R) into (random) direct orthogonal sum: 

n = ^nj, nr-=L\i^). (2.9) 



Correspondingly, let hoilf) denote the Friedrichs extension of operator to = —A/2 with do- 
main dom(to) = {f e L\I^) : f e €^(1^)}: 

{hn{I-)f){x):=-\{^f){x), (2.10) 

/ e dom(/i^(j;)) = {/ G W',{I^) : /(x7_,) = /(xj) = O} , 

where W2 denotes the corresponding Sobolev space. Then we get decompositions of the one- 
particle Luttinger-Sy Hamiltonian: 

hD = ®hD{I^), ujen, (2.11) 

j 

with domain 

dom(/i^) = domihoilf)) C H , (2.12) 



into random disjoint free Schrodinger operators [^ho{I'j) } . with Dirichlet boundary conditions 
at the end-points of the intervals {-^j^} .• The corresponding eigenfunctions have the form: 



^r,,D(^) = (0,0,...,^,^,^(x),0,...), (2.13) 

with eigenvalues {-E'sj(L^)} : 

h-^^>l^^ = E,^{L-)^>l^^. (2.14) 



Remark 2.2 For a given realization uo ^ Vt the spectrum of the random operator h2.10^ is 
explicitly defined by non- degenerate eigenvalues 

^(hnil-)) = [EsM^) = -^A , (2.15) 

with the corresponding eigenfunctions 

ijl^ix) = h.{x) MjM^ix - x-_,)) . (2.16) 

Here Ij'^lx) is the characteristic function of the interval I'^ . By consequence, the spectrum of 
the Luttinger-Sy Hamiltonian \2.11\) is the union of \2.1^] 

<^{hl) = \J<y{hD{I^)). (2.17) 

3 

By virtue of Remark 12.11 this spectrum is a.s. nonrandom, and it coincides with support 
of the integrated density of states A/". Moreover, in the case of the Luttinger-Sy Hamiltonian 
()2.1H) it is known explicitly since fTj. But the rigorous study and in particular the concept of 
" self-averaging" , which ensures this nonrandom property, are due to fH] • 



3 Self-averaging of the integrated density of states 

For the reader convenience we recall in this section some arguments that one uses to derive the 
spectral properties of the Luttinger-Sy one-particle Hamiltonian. Since our aim is to study the 
thermodynamic properties and Bose-Einstein condensation in this model, it is useful to derive 
the integrated density of states first for a finite system. 

3.1 Let A := [— L/2, L/2] C M. Then the finite Luttinger-Sy model with the Dirichlet bound- 
ary conditions at x = ±L/2 and with n — 1 singular point-repulsive (a — > -|-cx)) impurities 
corresponds to the one-dimensional self-adjoint Schrodinger operator 

n 

^L,x„:=0/iD(/,), (3.1) 

acting in the direct orthogonal sum of Hilbert spaces ()2.9p : 

n 

Ha:=0H,. (3.2) 

Here 

Xn = {xo = -L/2 < xi < X2 < ■ ■ ■ < Xn-i < x^ = L/2} , {Ij = (a;j_i,Xj)}"^^ , (3.3) 

and operators {/iD(-^j)}?=i are defined by ()2.10|) . 

To make this system disordered, Luttinger and Sy [TTj supposed that the impurity positions 
are random variables, which are independently and uniformly distributed over the interval A. 
Then instead of ()3.3p one gets the random sets {Xn}ujen '^^ (^' -^^ ^)' which a.s. contain n — 1 
ordered impurities {x'l-^ ._ . We denote the corresponding random Luttinger-Sy Hamiltonian 
and eigenf unctions in A by 

n 

hD,n,L ■■= hL,x, = hoilf) , h^^^^,^^;^^^ = E,^{L^)^^;^^^ , (3.4) 

where the eigenfunctions have the form: 

n 

^irAn = (O,O,...,^-,^(a:),...,O)G0H, , (3.5) 

see definitions fETT]l - (im|) . 

Notice that Remark 12.21 is valid in this case modulo the substitution of hf) by the random 
operator ()3.H) . In particular, for the spectrum of ()3.4p one gets representation: 

n 

The following proposition is an immediate consequence of the hypothesis about the inde- 
pendent uniform impurities distribution and the thermodynamic limit: L — > oo, n — > oo, with 
a fixed density of impurities 

A=lim ^. (3.7) 

L—fOO Li 



Proposition 3.1 (a) In the thermodynamic limit the above finite-volume random point field 
{X^} converges (in distribution) to the Poisson point field {X'^} with intensity A and corre- 
sponding random Poisson measure i2.^} . 

(b) The uniform and independent distribution of n — 1 points of impurities induces on A a 
random sets of intervals {if} ■_-. ,^ ^ ^, of random lengths {Lj = L^} . . The corresponding 
joint probability distribution is 

in - IV 
dPL^Li, -, L^) = \^J' 5{Lr + ... + L„ - L) dL^dL^ ...dL^. (3.8) 

(c) In the thermodynamic limit the lengths {-^^} . form an infinite set of independent random 
variables and distribution corresponding to \3.t^] converges (weakly) to the product-measure 
distribution a\ defined by the set of consistent marginals: 

k 
daxALn,. . . , LJ = X'H e-"^-rfL,, . (3.9) 

s=l 

The proof is standard 121221, see e.g. [121 ^^^ details. 

Recall that the finite-volume integrated density of states is defined by specific counting- 
function [201 • For operator ()3.4|) . it is a random variable of the form: 



1 ,^^ 1 /"^/^ 

^LiE):=j Yl 0{E-ElAn,L)) = - dx 9{E - hl^^^^){x,x) . (3.10) 



Here 6{E — h'^,^i^) (x, y) is the kernel of the spectral-projection operator oi h'^ ^ ^ corresponding 
to the half-line {—oo,E) and d{x) = I(o,+oo)(3;) stands for the step-function. 

Proposition 3.2 In thermodynamic limit the finite-volume integrated density of states i3.1U\} 
converges a.s. to non-random function 

with c = T\ I \pl. 

Proof: Explicit expressions ()2.15|1 and (j2.16|l imply for ()3.1()|1 the representation: 

A^r(.)4|:|.{.-(||)]. («.) 

Then by Proposition 13.11 and by (j3.7j) . (j3.12p we obtain 

Nm := a.^ lim K(J5) = a.s. lim - E E " ^ " (f ) = 

j = l S = l \ ^ ] ■' ) 

-^.{f<--(|)1|^A'Ef-.---(--(|)1. (-3) 



s=l \ ^ ' J ) s=\ 



The a.s. limit for (non-random) integrated density of states M\{E) exists by the Birkhoff 
ergodic theorem j2l 122] and the uniform convergence of the s - sum ensures the permutation of 
expectation with respect the ax - distribution (j3.9p and the sum. Thus, we obtain: 

Afx{E) = X^J2 rfL,e-^^' = A^e-^^^/^, (3.14) 

which yields the explicit formula ()3.11|) . D 

3.2 Formula ()3.1H1 allows us to recover for all energies E > the one-dimensional integrated 
density of states for the free operator t, i.e. the case when density of impurities A = 0, cf. 
Remark O and (jTTjl : 

/2 
limA6.(E) = A6.=o(^) = —Ve . (3.15) 

AjO TT 



Notice that for the Lebesgue- derivative n\{E) := dN'\{E)/dE, i.e. for the density of states 
Sect. 4, this limit is not uniform in E in the vicinity of the spectrum edge E = 0. This confirms 
the argument, previously presented in jTHj, that the Bose-Einstein condensation in such random 
media can not be viewed as a perturbation of the free case, since this phenomenon is tightly 
related to the behaviour of n\{E) near the edge [HIT^. 

On the other hand, for A > and for E close to the edge of the spectrum, the integrated 
density of states (j3.11|) exhibits the Lifshitz' tail behaviour: 

ATxiE) = \e-'^l^ {1 - 0(e-2^^/^)} , ^ i , (3.16) 

see Remark |2. 11 In this case Xvca^Ei^'^xiET) = 0. 

It is known ^TJ |TC] that behaviour (|3.16|) near the edge remains valid even if the parameter 
a > in ()2.6|) is finite. Notice that this parameter does not appear in the leading term of the 
asymptotics ()3.16|) . This can be explained by the fact that particle with small energy "sees" a 
point impurity potential with relative amplitude a/E ^ 1. Therefore, in spite of its singular 
nature the Luttinger-Sy Hamiltonian seems to be a good approximation for studying the BEG 
in Poisson random systems with non-singular repulsive impurities. 

4 Thermodynamics and Bose-Einstein Condensation 

The second quantization of the one-particle Luttinger-Sy Hamiltonian ()3.4|1 in the boson Fock 
space gives the one-dimensional PBG embedded into a random potential created by Poisson 
repulsive impurities ()2.6p with a = -(-cx). The latter implies that bosons are distributed over 
independent intervals ("boxes") {L'j^ . . 

4.1 We study the boson Luttinger-Sy model in the grand canonical ensemble, defined by the 
inverse temperature (3 > and the chemical potential /i. Since the model corresponds to 
independent "boxes" {L^} . , the grand partition function of the PBG in A = [— L/2,L/2] is 
the product of partition functions calculated in individual "boxes" : 



n 71 oo 

EtiP,^^) = l[El{P,f^) = l[l[{l-e 



-l3(E,ALf)-fi)' '^ 



see fl3.4|) . This gives for the grand canonical pressure 



p-(/?,;.) = -^5^5^1n(l-e-^(^».(-")-)) . (4.1) 

To ensure the convergence in ()4.H) we have to bound chemical potential from above: n < 
inis^^uj Es^{L'j). By virtue of ()2.15|) we obtain in the thermodynamic limit: 

a.s. lim iniEs^iL'^) =0 . (4.2) 

Lemma 4.1 For fi < and L ^ oo the pressure p%{j3,fi) converges almost surely to the 
non-random function 

p(/3, /i) = a.s. lim ptiP, /i) = --^ / dE nx{E) In (l - e''^^^-'^)) , (4.3) 

where the limiting density of states 

\2 pcX/VE 

nJE) := — . , (4.4) 

^ ' 2 ^3/2(ecA/^_l)2 ' ^ ) 

and c = n/V2, cf / fTTTl) . 

Proof: By definition of integrated density of states (j3.12p we can represent the pressure in (|4.1|) 
as the Lebesgue-Stieltjes integral 

pt{(3, /x) = -- J^ dUt{E) In (1 - e-^^^-'^)) . 

Then by virtue of (j3.13p we obtain that the limit 
p(/3,/i) = a.s. lim p'^(/3,/i) 

L—fOO 



\ 2 °° foo 

~ V / rfL,e-^^» In (l - e-/5((-/^»)^-^) 



exists a.s. and, after change of variables and calculation of the sum, takes the form: 

Setting k = y/E, we recover the relation ()4.3p with density of states ()4.4p . D 

Similarly we obtain the statement about the thermodynamic limit of the grand-canonical 
particle density. 

Lemma 4.2 For all fi < and /3 > 0, the grand- canonical particles density 

1 " °° 1 poo 1 

converges a.s. to 

nx{E) 
Jo 
with density of states n\{E) defined by (0^ 



P(/^'/^)= / dE ^^(^_;^ '_ ^ , (4.6) 
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Proof: By virtue of representation ()4.5|) . the demonstration follows the same line of reasoning 
as we used above for the limiting pressure (j4.3|) . D 

Corollary 4.1 By explicit formula ^-41 ) we obtain that for the Luttinger-Sy model, defined by 
the Hamiltonian \2.11\) . the critical density 

pjf3) = lim r dE-J^^4^ (4.7) 

is bounded. 

4.2 It is known that for PBG the condition pdP) < cxd implies the existence of (generalized) 
Rose condensation P, when the particles density p > pdP). To make it obvious in our case we 
have to study solutions p'^{P,p) of the finite-volume equations, see ()4.5|) : 

p = pt{P,p), cuen. (4.8) 

In fact, the asymptotic behaviour of /x'^(/3, p) studied for a general ergodic non-negative random 
potential in ^2]- These results then can be applied to the Luttinger-Sy model and lead to the 
following proposition: 

Proposition 4.1 Let p^{(3,p) be solution of the equation \4-^ for a given uj eVL. Then 

(a) for p < PciP) the limit 

a.s. lim /ii(/?,p) = pif3,p) < , (4.9) 

exists and is the unique root of equation defined by !i4.6\) : 

p = p{(3,p), (4.10) 

(b) for p > Pc{(3) the limit 

a.s. lim /ii(/5,p) = , (4.11) 

L^oo 

For p > PciP) the limit (jLlip implies that the density of condensed particles can be define 
in the usual (for generalized condensation) way: 

Po(/3,p) :=lim |a... \un^ J^^idE) ^^^^_^^^^^^^^^ _ ^ ^ = p - p^P) , (4.12) 

see e.g. 0. If p < Pc{/3), the limit is zero. Notice that this nonrandom limit is a consequence 
of the chemical potential asymptotics (Proposition 14. ip and of the uniform convergence of the 
particle density ()4.5p . see [12], Theorem 4.1. 

5 OfF-Diagonal Long-Range Order 

In this section we study the problem the two-point correlation function [^ |21] . By definition, 
in the finite volume A and for any cj e f2, it has the form: 

n oo ^ 

j = l Sj=l 

n oo _ 

= EE .(.,.(.p-.)_. ^(^)€.(^)^ (5-1) 

j = l Sj=l 

10 



where (■ , ■)Rn denotes the scalar product in W^, see fj3.4|) and ()3.5p . Therefore ()5.1|) is the kernel 
of for the one-body reduced density matrix, see e.g. |15j . 

5.1 We know that this function is not self- averaging in the thermodynamic limit P^ I2UJ. To 
get a way out, we proposed in jTS| to consider the space-averaged version of (|5.H) : 



1 /"^/^ 
Pu{x,y;P,fj.) := - dz p^{P, fx;x + z,y + z) . (5.2) 



The motivation was based on the fact that in the limit A | 0, we recover from (|5.2j) the free 
case, see jTSj and Section 5.2 below. 

Theorem 5.1 For the Luttinger-Sy model, the thermodynamic limit of the space- averaged two- 
point correlation function i5.^} for (3 > and fi < 0, is a.s. nonrandom and has the form: 

p{x, y; /3, /x) = po(/5, p) + e-'\^-y\ j^ dE-^^^ cos{V2E{x - y)) . (5.3) 

Here nx{E) is defined by \4--4\l (ind po{P,p) the condensate density \4-1^4l - 

Proof: We consider first the case of negative chemical potential ()4.9|1 . i.e. po{P,p) = 0. Using 
explicit form of eigenf unctions ()2.161 we obtain for the thermodynamic limit of ()5.2j) : 

, n oo ^ 

p(x, y; /?, p) = ^im - E E ^/.(i..,(Ly)-,) _ ^ >< (5-4) 



j=l s=l 



2 '•^ 



rfasin(/c^^,(Lp(x + a - y^'^i)) sin(/c^^.(LJ)(y + a - i/^_i))I/-(x + a)Ij^{y + a) 
-^i Jo 



with ks^iL'j) := ./2Esj{L'j). Let us put, for simplicity, x > y. Then the integration is reduced 

to [yj — y,y'^ — X -\- L^], such that (x — y) < L^. Since ks^{L'^)L'^ = stt, the integration over a 
yields 



|- /" dasm{ks^{L'^){x + a- yJ^J) sin(A;,^.(L5^)(y + a - y}^_i))%(a; + a)I,^(y + a) 
cos{h^{L-){x - y))e{L- -{x- y)) (l - ^ 

with step function 9{z). Since by Proposition 13. II random variables L^ are independent in the 
limit L ^ oo, we apply to ()5.4j) the Birkhoff ergodic theorem and find the limit: 



p{x,y; (3, p) = a.s. A^ >^ / dLjC' 

„_i "'0 



cos(y2cs(x - y)/Lj)9{Lj -(x-y)) [l 
with c = 7r/v2. If we put q = cs/Lj, then 



' g/3((cVi,)2)-M) _ 1 ^ 



J' 

oo ^oo 



/"OO 7 

p{x,y;P,p) = A^ V / — 



-csX/q 



e/^iq'~i^) _ 1 



X (5.5) 



cos(V2g(x - y)) 6{s - q{x - y)/c) c {s - q{x - y)/c} . 
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The sum over s yields 

~ -cA/g 

E '^"'"''('^ - ^(^ - yy^) = '-''^-' \i_,-..„Y ' (5-6) 

where Smin = [(lix — y)/c] denotes the entire part of q{x — y)/c. Then after change of variables, 
q = y/E, we find by (jOD and (jS^D for /i < 0: 

p{x, y- /3, /x) = e-^l^-^l T rfE ^('^"Jj^|_ ^ cos(v^(x - y)) . (5.7) 

We put here \x — y\, since the proof for x — y < is identical to that for Q < x — y. 

Now we shall study the case when the condensate exists. The finite-volume chemical po- 
tential iif^{P,p) is a solution of equation (j4.8p . By Proposition I4.ir b) for p > pdP) it implies 
(I4.11|) . i.e. p'^{P,p > Pc{P)) converges a.s. to 0. To find the limit of space-averaged correlation 
function ()5.2|) for the sequence {pLiP^P)}L ^^ ^P^i^ ()5.2p into two parts: 



-. n oo -, 

j=l s=l 

COS [,Je^){x - y)) e{L- -{x- y)) {l - ^^) e{5 - E,^{L-)) 

_ n oo 



j=l s=l 



COS [.Je^){x - y)) e{q - (x - y)) {l - ^^) e{E,^{q) - 5) , 

for some 5 > Q. Since in the second term of the right-hand side of (j5.8p the eigenvalues 
Es-{L'^) > 5 > 0, the limit ()4.1ip and uniform convergence of the sums (cf. Corollary 3.1 in 
di) yields 



X 



-. V CO ^ 

lim lim pfiix^y: 3, piiB, p)) = a.s. lim lim — > > ^,^ ,,^. — zrnr-V\ 

j=l s=l 

COS [,/e^){x - y)) em -{x- y)) {l - ^^^) e{6 - E,^{L-)) 

+ p(x,|/;/3,0), (5.9) 

where the last term is defined by ()5.7p . 

To study the limit of the first term in the right-hand side of ()5.9|) we use the fact that the 
levels with energies Esj{L'j) < 6 are, by definition, of order 0{6). By p.l5|l these levels are 
defined in large "boxes" of lengths of order 0{S~^^'^). Then for Es.{W) < 5 and |a; — 1/| -C 5~^/^ 
with S small enough, we obtain asymptotics 

cos [^/e^){x - y)) e{L- -{x- y)) {l - ^^^) = 1 + O (v^) . 

Therefore, by definition of generalized condensation ()4.12|1 and by (j5.9|l . we get 

p^(x, y- 13, 0) := a.s. lim lim p^(x, y- (3, pl{l3, p)) = po(/3, p) + p(x, y] (3, 0) , (5.10) 

5J.0 L^oo 
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for p > PciP)- This finishes the proof of ()5.3|) . D 

5.2 It is instructive to make a contact between the concept of the space- averaged one-body 
reduced density matrix (j5.2j) and the one for nonrandom free PBG. 

Corollary 5.1 When density of particles p exceeds the critical value pdP) d^.Tp , the space- 
averaged one-body reduced density matrix of the Luttinger-Sy model, manifests Off-Diagonal 
Long-Range Order: 

0DLR0{f3,p):= lim p(x,|/; A/i(/?,p)) = po(/5,p) , (5.11) 

\x-y\^co 

We see that similarly to the nonrandom case this limit is defined by the condensation density 

urm . 

Remark 5.1 Notice that for p < (or p < Pc{/3)) , the space-averaged reduced density matrix 
/ I5.5J) remains consistent with the free nonrandom case. Indeed, by ^■4\ ) we get for (|5. 7[ ) the 
limit 

\imp{x,y;/3,p) = p{x,y;f3,p) , (5.12) 

AJ,0 

where 

p{x, y- P,p) = - J^ -^ g^(i;-^) _ ^ cos(y2^(x - y)) (5.13) 

coincides with two-point correlation function of the free PBG 11^ . This equivalence is valid 
only when there is no condensation, since by ^-4] ) o-nd (|^.7| ) one has that lim^j^o Pc(/^) = oo. 

Let the single-impurity potential u{x) > 0, a; G M, have compact support [—77/2,77/2], see 
Section 2.2. Then one has 

Proposition 5.1 [113] Let 7 := 1 — e""", with v := j^v{x)dx. Then 

p(x, y- /3, p) < p(x, y- /?, p) e-^^d^-^l"") , (5.14) 

P— a.s. for any p < 0. 

Since for the Luttinger-Sy model the single-impurity potential is defined as a singular (5-point 
potential with infinite amplitude (Section 2.4), we get in ()5.14j) 7 = 1 and 77 = 0. 

To check this directly, put E = w\x — yl"^. Then for p < we can represent two-point 
correlation function (|5.7p as 

'^ POO 

p(x, y; (3, p) = Re e"^!"-^! V e^''' / rfEe-^^^r2A(^)e'^^l"-^l (5.15) 

= e-^l"-^l \x - y|2 Re Y^ e^^' / dw nx{w\x - yp) e-l-s^l'(^-^-^^^) . 

.=1 J^ 

To calculate the asymptotics of p{x,y; f], p), when |x — ?/| — >■ 00, we estimate the last integral 
in ()5.15|) by the saddle-point method. Then ()4.4j) implies 



°° g-|x-y|2/4/3s 



p{x,y;P,p) = e-^l-^l O^e^^'^— -— -^ + e-V%ll-^IO(|x-y 



_1_ 1^ V ■■■■ "'1/11/ II \ 

(27r/5s)i/2 



=1 
-^l^-^l p{x,y;P,p) + e-(v^+^)l--^l 0{\x-yn , 



for \x — y\ -^ 00, with the free PBG two-point correlation function p{x,y] (3, p) defined by 
fl5.13j) . This confirms the statement ()5.14j) for the case of the Luttinger-Sy model. 
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6 Comraents and discussion 

6.1 Critical density 

We start by remark concerning modifications of tlie Luttinger-Sy model properties, and in 
particular of the value of the critical density, when one passes from infinite to a finite amplitude 
a < oo of the ^-potential ()2.6p . 

Recall that operators {h^}g^^Q correspond to a monotonically increasing family of quadratic 
forms with /ia=+oo = ^z)' ^ee ()2.^|) . Then by definition of the integrated density of states ()3.1U|) . 
(j3.11|) and by the mini-max principle for U^ and U^^ one gets 



Nx{E) = Afx,a=+o.{E) < f/xAE) < Arx,a=o{E) = Afx=oAE) = ^ ■ (6.1) 

Notice that the integrated density of states for the free case a = coincides with that for the 
zero impurity density ()3.15|) . A = 0. From ()6.ip one gets the corresponding inequalities for 
critical densities ()4.7|1 indicating the enhancement of the BEC, that was remarked already by 
Kac and Luttiger |H1 E] • 

More refined arguments (see e.g. [20], Ch.III, 6B) give for E < 7r^a^/32, the estimate, cf. 
(HH: 

oo 
Mx{E) < Ux,a{E) < A ^ e-A(c/v/S-4/a) ^. ^^*^^^) _ (g_2) 

Then, definition of the critical density ()4.7|) and ()6.H) yield 



"'l^'^)^ = Sfgls^^f ''^^^W 



e/3E _ 1)2 

< j dE AfxAE) .Jj, _ ^.2 =■■ PcW, A, a) , (6.3) 

with obvious limits: lim.a^o Pc{P , A, a) = oo and limx^o pdP , A, a) = oo, by (j6.ip . 

The estimate (j6.2p yields the upper bound on Pc{/3, A, a) for small a. Setting E{a) := (7ra/8)^ 
by (jHUll-dniSl) we get that 



Pci^^ ^' «) < I dE ATl^iE) ^^^3^ + y^ ^ dE^ ^Jl _ ^^^ =: J(/3, A, a) . (6.4) 



lEia) 

Then for fixed A > and small a > we obtain the estimate: 



p,(/3,A,a)</(/5,A,a)<-i-f-) 7-7^ + -^- 

pA Vtt/ 4e(v2 - 1) a (Jir^ 

Similarly, for a > and small A > we obtain by ()6.4j) that 

,^ , , 1 /-^("^ dx e-^/v^ r ,r. v^ Pe^^ 

^^^^' '' ^^ - A^ I ^T^^^ + L ^ V^^^ • 

Notice that the bounded critical density pdP, A, a) for A > and a > is the key criterium 
of existence of BEC in the one-dimensional system with integrated density of states Mx,a{E), 
cf. ()4.12j) . On the other hand the Bogoliubov-Hohenberg theorem says that there is no BEC in 
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translation invariant boson systems if dimension is less than d = 2, see e.g. 013]. Therefore, 
the BEC in the Luttinger-Sy model is of a different nature than the case without random 
impurity potential. 

In fact, the randomness of the impurity potential is not indispensable for BEC in one- 
dimensional perfect Bose-gas. To this end we construct nonrandom hierarchical models with 
impurity potential which manifest the BEC via mechanism similar to that in the Luttinger-Sy 
model. 

6.2 Hierarchical model for BEC in one-dimensional nonrandom intervals 

We present here a nonrandom hierarchical one-dimensional system, which manifests BEC and 
in a certain sense mimics the Luttinger-Sy model. 

Type I BEC. Let A := (0, L) be a segment separated into n impenetrable intervals of lengths 
Lj, j = 1, ..., n such that A = n/L < cxd. For simplicity we take the hierarchy when all intervals, 
except the first {largest) one, are identical: 

Li = — - — and Lj^i = Ln = — . (6.5) 

A n — 1 

Then one gets 

lim Li = +00 and lim Ln = -r ■ (6.6) 

L— >oo L— >oo A 

This non-random system presents an obvious analogue of the Luttinger-Sy model. Here again, 
the quantum states are defined in independent intervals and they have energies 

^hs = -pr ' 3 = ^,--,n , s = 1,2,... , (6.7) 

with c^ = 7r^/2. The spectrum of the corresponding Schrodinger operator is discrete and 
bounded below by zero, cf. (j2.15|l . ()2.17j) . Then the chemical potential is // < and the PBG 
particle density in A has the same expression as in ()4.5|) 

PL(A/i) = iY.T. em.s-,) _ ^ . /^ > . /^ < • 

j=l s=l 

By virtue of the hierarchical structure of intervals we can separate the expression for density 
into two parts: 

1°° 1 n-l°° 1 

s=l s=l 

Since Li = 0(ln(AL)), the first sum in ()6.8|) converges, when L — ;> oo, to zero for all yU < 0, i.e. 
we obtain 

oo ^ 

p(/3,/i) = ^im pi(/3,/i) = >^J2 ^mcs/x)^',) - I ■ ^6-9) 

As a consequence, the critical density for this system is finite: 

oo _ 

PciP) := sup p(/5, /i) = p(/3, 0) = A ^ ^pics/xr - 1 ^ °° ' ^^-^^"^ 

^— s=l 
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and we have BEC condensation, when p > pdP). 

This condensation is of type I, since the difference between the ground-state energy and the 
energy of the first excited state (which are both locahzed in the biggest interval Li) is of the 
order 0(L|f^) = 0((ln(AL))"^), see e.g. [T], or j23]. In this case the solution pl(P,p) of the 
equation 



P 



1 1 >^ 1 



S>1 



V — 1 sr^ 

+ ^E 



L ^^ g/3((cs/L„)'-M) 
s=l 



for p > Pc{f3) and large L has asymptotics 



Pl{P, p) = i?i,i - . ,1 .. + 0(1/1') , (6.12) 

PPo{p,P)L 



see (16. 7p . Inserting (J6.12J) into (J6.1H) we obtain in the limit 



^ 1 ^ l_ 



= Po(/5,p) + Pc(/5) 
where by fl6.12|) the condensate density is 



Po(/5,p) = lim - 



L^oo L e/3(cVi?-ML(/3,p)) _ 1 ■ 

Therefore, this one-dimensional hierarchical model shows a type I BEC localized in the [logarith- 
mically) large, but not macroscopic, domain corresponding to the ground-state wave function. 
Generalizations to another hierarchy of intervals with one largest interval trapping BEC are 
obvious. 

For example, it is easy to generalize the above observation to the type I BEC localized in a 
finite number of M identical (logarithmically) large intervals: 

L, = ... = Lm = -^^ and Lj = L„ = -/■ , M<j<n, 6.13 

A n — M 

cf. (16. 5j) . Then similar to the case M = 1 (j6.8|) one gets 

1 Y^ Y^ 1 n — M Y^ 1 

which implies trough verbatim that for M > 1 the critical density (J6.1U|) rests the same. If 
p > PciP), the equation p = pl{j3,p) (j6.14p yields for asymptotics of the solution pl{P,p) an 
expression similar to ()6.12|) for M = 1: 

'"('^•'') = ^» - ft;^ + «(^/^'> . 1 < . < M . (6.15) 
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Here £'1,1 = . . . = Em,i by dHUI) and ^J^ . Then taking into account dHUj) and ^J^ . (IFTT^ 
we obtain by (j6.14|) 



M 
i ' 

P 






= Po(/^,p)+Pc(/?) , 
where the condensate density is equally shared among the first M intervals: 

M 1 

Po{P,p) = lim — 



Type II BEC. To obtain the type II BEC in one interval we take, instead of ()6.5|) : 

Li := a/L/A and L^yi = L„ = — . (6.16) 

Then ()6.8p gets the form 

Pl(/3, p) = ^ E e/3(Ac^.VL-M) _ 1 + ^ E ,,i^syLi-^ _ 1 • (6.17) 



s=l 



Since for fi < the first sum in ()6.17|1 converges to zero, when L -^ 00, we obtain for p(/5, /i) = 
liniL^^oo Pl(/5, /^) and pdP) the same expressions as in ()6.9p and ()6.1()j) . 

Now, if p > Pc{P), the solution iil{(3,p) of equation p = pl{P,p) has asymptotics defined 
by dnHI): 

-.00 ^ 00 -. 

s=l s=l 

= po(/3,p)+p.(/3). (6.18) 



As in (I6.12J) this implies 



PlW, p) = i^i.i - ^^^ + 0{l/L') , (6.19) 



see fj6.7p . where by ()6.18|) the coefficient A{[3,p) > satisfies the equation 

00 ^ 

Hence, for p > pc(/?) the BEC 

00 -. 

is splitted between infinitely many states in the largest interval Li, i.e. this is condensation of 
the type II, p. 
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Type III BEC. Now we show that (unusual) spatially fragmented type III BEC is possible in 
our hierarchical model. To split BEC between infinitely many states in different intervals, let 
volume A be occupied by [ln(?7, + 1)] identical [logarithmically) large intervals: 

^_HXL)_ ^ l<J<Nn + l)]=:M„ and L,>m„ = ^n := ^J"_^^" . (6-21) 

for Mn < j < n, of small intervals, cf. (jH^l)- Then (similar to ()6.14j) ) we get for the particle 
density 

1 ^-^ Y^ 1 n — Mn \r-^ 1 

PLiP, /^) = ^ L L ^.ic^syL^,) _ 1 + -^— L ,,i.^syLi^,) _ 1 • (6.22) 

j=l s=l s=l 

Since by (I6.2ip we have limi^oo -^n = liniL_^oo (^ — Mn)/ L = A, then the critical density (|6.1(J|) 
rests the same. If p > pdP), the equation p = pl{[3,ij) (|6.22|) yields for asymptotics of the 
solution pl{P, p): 

Then we obtain for the particle density in large intervals 

1^ ^ 1 ^ Mn 1 

i = l 5 = 1^^ i 

M °° 1 

+ il^T^E e/^((Acs/lnnP-M,(/3,p))_l = ^o(Ap). (6.24) 

s=2 

The limit ()6.24|1 and ()6.22|1 . ()6.23|1 imply that the condensate density po{P,p) = p — pdP) is 
splitted between ground states of infinitely many logarithmic intervals in such a way that the 
condensate density in each interval is zero. This is a spatially fragmented type III BEC, which 
is different from the that corresponding to spread out over infinite number of states discussed, 
e.g., in [11123]. 
6.3 Statistics of large Poisson intervals 

By virtue of 6.2 to discriminate between possible types of BEC in the Luttinger-Sy model 
one has to study statistics of the size {L'j} ■ of intervals induced by Poisson distributed point 
impurities, see Proposition 13.11 

In fact, the first attempt to elucidated this question is already contained in ^H]- They gave 
some arguments in favour of that for large finite A the largest interval Ii has a typical length 
of the logarithmic order: 

L^ ~ A~^ In(AL) , L ^ oo . (6.25) 

Moreover, along the same line of reasoning they conclude that all other intervals I'^^-^ are 
typically much smaller than I^ . Then neglecting the fluctuations the length of L^ they conclude 
that BEC has to follow scenario we described in 6.2 as type I BEC, cf. ()6.5p and ()6.25|) . 

To check these arguments and to bolster them by some rigorous reasonings we use Propo- 
sition 13.11 First we note that the average length of the Poisson intervals is 

POO 

E^^{Ll) = \ dLLe^^^ = \^^ , (6.26) 

Jo 
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i.e., the total average length of any sample of intervals {if} ._, is k/\. We are interested into 



3 i i=i 



> 



density of the joint probability distributions generated by the events: {u E Vt : Ll-^ > L 
... > L'^ }■ They have evidently the form 

d(^xALj„ ..., LjJ := k\ 9{Lj^ - Lj^) 9{Lj^ - Lj^)...9{Lj^_^ - Lj^) d(rx,k{Ln, • • • , LjJ . (6.27) 

Then one gets for the joint probability density of the two largest intervals: 

dal,,{Lj^,Lj^)/dLj^dLj^ = k{k - 1) A^ e-^^^^e'^^^^ (l - e-^^^^)'"^ e{Lj^ - L^J . (6.28) 
Now, let A{s,t) be defined for non-negative integers s and t by 

A{s,t):= dxln'{x)x^e-'' . (6.29) 

Jo 

For s > 1 and t > 1 this function verifies the following identities: 

A{s, t) = tA{s, t-l) + sA{s - 1, t - 1) , A{1, t) = tA{l, t - 1) + T{t) , (6.30) 

where r(t) stands for the Gamma- function. Then by ()6.28p and ()6.29p . ()6.30|) we obtain for 
expectations of the two largest intervals: 






s=l 

k 



(-«)-4^-^4Ei. (-) 



"^.*' ^^' k\\ A X^s 



s=2 



By virtue of (16.311) and (16.321) the mean difference Eg.> (L^ — L^) = 1/A is independent of the 

X,k -^l -J ■^ 

number k of intervals in the sample, whereas they have, for large k, the logarithmic size (cf. 
(E2SD): 

K>JLI^^) = J Hk) + ipi,2 + Oil/k) , (6.33) 

with respect the total average sample length k/X, here Pi = C := 0,577..., is the Euler 
constant, and P2 = C — 1. Using ()6.27p and ()6.29p . ()6.30|) we find that the variance of the 
difference between two largest intervals in the sample is also k-independent and has the form: 

^«V.(^i.-^T.) = ^- (6.34) 

Moreover, by the joint probability distribution ()6.28|) we obtain for any S > that probability 

P{c^ : LI - LI >6}= e-^^ (6.35) 



of the events Ak{6) = {u : L^^ ~^'j2 > ^}) is independent of k for increasing sequence of samples 
\lf\ . ,, when k ^ 00. 

By 6.2 and ()6.33p . ()6.34p we see that the type II or /// BEG are impossible in the one largest 
logarithmic "box", since that total average length of the sample is k/X. To exclude the type 
I, II, III condensations via a space fragmentation between, e.g., two "boxes" , we have to estimate 
probability of events corresponding to the state-energy spacings between two largest intervals. 
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By 6.2 (see fl6.15|) . ()6.19p . ()6.23|) and J2]) this spacing should be larger than inverse of the 
total sample length, which is {k/X)~^. To this end it is sufficient to estimate the probability of 
the event Sk{a > 0,7 > 0) corresponding the spacing between ground states: 

P{5,(a,7)} :=P{^ : i?.=i(L- (A;)) - E.=i(L- (fc)) > ^} . (6.36) 

Here we denote the energies in the sample [if] ._-. by 

EsiLlik)) = j^^^ , r = l,...,k, s = 1,2,... . (6.37) 

Notice that by ()6.35|) we obtain that there is a kind "repulsion" between energy levels in different 
intervals. Indeed, 

P{^,(a,7)} > n^ : Ll{k)-Ll{k) > ,^^^Ll{k){Ll{k)f} = 



00 /"OO 



where limfc^ooPfc(fl, < 7 < 1) = 1 by explicit calculations in ()6.38p . The same argument is 
valid for other than ground states as well as for intervals ilf\ „ instead of If . Therefore, 
in this limit with the probability 1 the spacing is too large for fragmentation of condensate 
between the largest and other intervals. 
6.4 The Kac-Luttinger conjecture 

The above arguments prove the Kac-Luttinger conjecture in the case of the one-dimensional 
random Poisson potential of point impurities: for PBG the BEG is of type I and it is localized 
in one "largest box". 

To make this statement more precise recall that BEG exists only in the thermodynamic 
limit, which we construct as increasing family of samples of intervals {if} ._, induced by the 
point impurities on M. Since these random variables are independent, we can choose increasing 
sequence of independent samples with one largest interval and with the property that 

k 

lim VP{^fc(a,7)} = 00 , (6.39) 

r=l 

see f)6.38|) . Then by the Borel-Gantelli lemma j22] 

P{lES^fc(a,7)} = l, (6.40) 

where the event 

00 

lim 5*^(0,7) := p] y 5/(0,7) 

k=ll=k 

means that infinitely many events {Sk{a,'y)}k>i take place. Together with 6.3 the statement 
(J6.4UJ) mean that with probability 1 in the thermodynamic limit M the BEG is localized in a 
single "largest box", and this condensation is of the type I. 
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